Silicon Valley
tydiverse
tydi_tools
utopia
sf
dib <- read.table(
file = "data/dib2020.txt",
header = TRUE,
na.strings = "NA"
)
summary(dib$Kda)
## Min. 1st Qu. Median Mean 3rd Qu. Max. NA's
## 0.370 1.445 2.650 2.620 3.485 11.230 199
library(tmap)
library(sf)
## Linking to GEOS 3.8.0, GDAL 3.0.4, PROJ 6.3.1
library(tmaptools)
tmap_mode("view")
## tmap mode set to interactive viewing
dib_sf <- st_as_sf(dib, coords=c("X.UTM20","Y.UTM20"), crs=32720)
tm_shape(dib_sf,
name = "Unidad cartográfica") +
tm_dots(
col = "Kda",
title = "Kd atrazina",
palette = "-viridis",
popup.format = list(
digits = 1,
decimal.mark = ",",
big.mark = ".",
suffix = " L/Kg"
)
)
littleyoda
Material en: